Electronic Heat Transport Across a Molecular Wire: 
Power Spectrum of Heat Fluctuations 
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With this study we analyze the fluctuations of an electronic lieat current across a molecular wire. 
The wire is composed of a single energy level which connects two leads which are held at different 
temperatures. By use of the Green function method we derive an explicit expression for the power 
spectral density of the emerging heat noise. This result assumes a form that is quite distinct from the 
power spectral density of the accompanying electric current noise. The complex expression simplifies 
considerably in the limit of zero frequency, yielding the heat noise intensity. The heat noise spectral 
density still depends on the frequency in the zero-temperature limit, assuming different asymptotic 
behaviors in the low- and high-frequency regions. These findings evidence that heat transport across 
molecular junctions can exhibit a rich structure beyond the common behavior which emerges in the 
linear response limit. 

PACS numbers: 05.60.Gg, 73.63.-b, 68.65.-k 
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I. INTRODUCTION 

The experimental realization of molecular junctions [l| 
some fifteen years ago has sparked a wave of research ac- 
tivities both in theoretical 0t3 ^-^id the experimental 
[3-[g] communities. Presently, single molecule electron- 
ics is considered as a possible potential substitution of 
the silicon-based elements in the information processing 
technology 0, Q , and this perspective makes studies of 
molecular transport properties very appealing. Besides 
the standard current-voltage characteristics [l|, 0, Q, it 
is also possible, for example by use of the full counting 
statistics [9l-[ll| , to obtain the information about the fluc- 
tuations of the electric current flowing through a molec- 
ular wire dMl. 

The issue of heat transport across such molecular junc- 
tions attracted much less attention in the prior litera- 
ture. This may be so because a measurement of heat flow 
experimentally is far from being a straightforward task. 
The problem of heat current constitutes, however, an im- 
portant physical issue. This is so because the structural 
stability of any molecular structure depends sensitively 
on the heat flow which accompanies the inter-electrode 
charge transport. 

With the systems of interest being of nanoscale the 
corresponding heat noise can also be quite large. This 
may be so even in situations where the average heat cur- 
rent is identically vanishing all together, as it is the case 
in thermal equilibrium. Moreover, the properties of noise 
correlation features, or likewise, its frequency-dependent 
spectral properties and, as well, its zero-frequency power 
spectrum, are in no way directly related to the mean 
value of the heat flow itself. It is thus of outmost interest 
to gain some insight into the size of heat noise in such 
molecular wire setups. In particular, it would be very 
useful to have analytical estimates for its power spectral 
density (PSD) available, even at the expense that these 
may predominantly apply to idealized setups only. 

Heat transport across conducting molecular wires, ex- 



plicitly induced by a difference of the lead temperatures, 
Tl and Tr, is a phenomenon which involves electrons, 
phonons and their mutual interaction [J, [17h25| . |27| . 
Therefore, the definition of heat, carried through the 
wire, should be addressed with care, with the need to 
distinguish between heat transfer mediated by electrons 
and a one mediated by phonons. 

The latter issue constitutes the realm oi phononics [28| . 
a promising novel research area which may lead to new el- 
ements, such as molecular thermal diodes, thermal tran- 
sistors, thermal logic gates, to name but a few [2^, [28l - 
|3]| . Here also, the size of fluctuations in heat current 
does matter; this is so because those may well turn out to 
be deleterious to intended information processing tasks. 

Energy transport mediated by electrons is a process 
which is related to its electric current: electrons are mov- 
ing from lead to lead, carrying not only charge but also 
energy. However, the amount of energy carried by a sin- 
gle electron, unlike to its charge, is not constant [3j, [SJ] . 
In contrast to the studies that examine the average heat 
flow much less attention, however, has been paid to the 
issue and impact of fluctuations of heat flow across vari- 
ous nanoscale devices. In some prior work |35| the heat 
transport through a ballistic quantum wire has been con- 
sidered in the Luttinger-liquid limit, by neglecting the 
discreteness of the wire's energy spectrum. In more re- 
cent publications, Rcfs. [Sa, |33], the PSD of the heat 
current fluctuations has been derived within the scatter- 
ing theory, under the assumption that the electrons are 
transmitted (reflected) at the same rates, independently 
of their actual energies. The results of the last two pa- 
pers are distinctive because they have shown that the 
heat noise exhibits a well-pronounced frequency depen- 
dence even in the zcro-tcmpcrature limit. 

With this work we consider the electronic heat cur- 
rent that proceeds across a molecular wire composed of 
a single energy level with the two leads held at a con- 
stant temperature difference. In doing so we shall ne- 
glect electron-phonon interactions and electron-electron 



interactions. Such a simplification can be justified in sit- 
uations that involve a short wire. Then, the Coulomb 
interaction via a double occupancy shifts the energy far 
above the Fermi level so that its role for thermal trans- 
port can be neglected. Likewise, the electron dwell time 
is very short as compared to the electron-phonon relax- 
ation time scale. 

In contrast to prior works [3 a . |37[, however, we here 
take into account the dependence of the transmission co- 
efficient on its electron energies, and, within the nonequi- 
librium Green function approach 0, [2^, derive an ex- 
plicit expression for the PSD of the heat current fluctua- 
tions, S^{u}). With this result at hand we explore differ- 
ent regimes of electron transport and demonstrate that 
the heat noise in fact is quite distinct from its electric 
counterpart. 

In Sec. Ill, we evaluate by use of the Green function 
method the expression for both, electronic heat current 
and its noise. We study analytically the distinct proper- 
ties of the electronic heat noise and compare this result 
with the electronic current noise counterpart. In doing 
so we demonstrate that for the chosen setup of a molecu- 
lar junction with a single molecular energy level different 
power laws for the PSD occur in different frequency re- 
gions. We discuss the obtained findings and end with 
conclusions in Sec. IV. 



II. MODEL SETUP 

The molecular junction setup is depicted with Fig. [T] 
It is described by a Hamiltonian 



H — H^ 



H, 



leads 



He, 



(1) 



containing three different contributions, namely the wire 
Hamiltonian, the role of leads and the wire-lead coupling, 
respectively. We consider here the regime of coherent 
quantum transport whereby neglecting dissipation inside 
the wire. The wire is composed of a single orbital; i.e.. 



H^ 



(2) 




at an energy Eq, with the fermionic creation and an- 
nihilation operators, d^ and d. The energy level Eq 
can be tuned by applying a gate voltage. This ide- 
alized setup makes possible explicit analytical calcula- 
tions. It mimics a double barrier resonant tunneling 
structure GaAs/Al^Gai -^-structure of the type consid- 
ered for electronic shot noise calculations in Ref. |3g |. 
herein truncated to a single Landau level. The leads are 
conventionally modeled by reservoirs, composed of ideal 
electron gases, i.e.. 



H, 



leads — / ^ ^^Q^lq^^Q 



(3) 



where the operator cl{c(q) creates (annihilates) an elec- 
tron with momentum q in the (. =L (left) or £ =R (right) 
lead. We assume that the electron distributions in the 
leads are described by the grand canonical ensembles 
at the temperatures Tl/r, and with chemical potentials 
/iL/Pj. With such ideal electron reservoirs we then obtain 

(cJ^Q'g') = Su>6qq'fi{eig), where 



fd^lq) 



o(siq-fJ,e)/kBTi I 2 



(4) 



denotes the Fermi function. 

Next we impose a finite temperature difference AT = 
Tl — Tr and use identical chemical potentials, Hh = ^J'R — 
H for the leads. When an electron tunnels out from a lead, 
the energy E leaks into the wire. This energy presents 
a heat contribution, 5Q, which in terms of the chemical 
potential, /i, reads 5Q = E — ^i. In the following we 
assume that all the electron energies are counted from 
the chemical potential value /i, i. e., we set ^ = |4l.[l7|. 

The tunnelling Hamiltonian, 



He, 



tq 



d + h.c. 



(5) 



mediates the coupling between the wire and the leads. 
The notation h.c. denotes Hermitian conjugate. The 
quantity Ve^q is the tunnelling matrix element, and the 
tunneling coupling is characterized in general by a spec- 
tral density, Tt{E) = I'nY.q l^«9p'5(^ - ^(-q) II- In what 
follows we also use the wide-band limit of the lead con- 
duction bands, with Y i{E) = F^ . 



III. FLUCTUATIONS OF HEAT FLOW 



FIG. 1: (color online) Setup of a molecular junction: two 
metal leads, filled with electron gas, are connected by a sin- 
gle orbital eo- The coupling strengths are determined by the 
constants Fl/h. The left lead is prepared at a higher temper- 
ature as compared to the opposite right lead, i.e. Tl > Tr,. 
The chemical potential, \i,, is the same for both leads so that 
no electric current due to a voltage bias is present. 



We perform the derivation of the central quantities, i.e. 
the heat current and its PSD of corresponding fluctua- 
tions, within the Heisenberg description. Setting for the 
energy operator. 



Ex 



— 2_^ ^LqCi^qCl^q , 



(6) 



its time derivative yields the operator for the heat cur- 
rent, reading: 



2e 



Jtit) = -Y.^iHVi.A,(mt)] 



(7) 



This current is positive when heat transport proceeds 
from the left (i.e. hot) lead to the right (cold) lead, see 
in Fig. 1. Upon formally solving the Heisenberg equation 
of motion for the lead operators, we obtain 



aS) =c,,(to)e-*-''« (*-*»)/" 



iVlq 



^ / di'e-*^''(*-*')/''d(i') 



(8) 



where the first part on the right hand side describes the 
dynamics of the free electrons in the leads, while the sec- 
ond part accounts for the influence of the wire. 

By inserting Eq. ^ into the Heisenberg equation for 
the electron annihilation operator within the wire, we 
find 

rf= j£orf-^^i-^d + a(t)+CRW, (9) 



where 



-T^lq{t - to) 



Clgito) . (10) 



To obtain the solution of Eq. ©, we follow the Green 
function approach in Ref . [j] , and first solve the follow- 
ing differential equation 



(! + !££ 



2h 



^)Git-t') = Sit-t'), (11) 



and then apply the convolution d{t) = / G{t — t'){£^i^{t') + 
^R{t'))dt'. The solution of Eq. ([TT]) is given by: 

G{t) = 6l(t)e-"^«*/''-(^L+^«^*/2'* . (12) 

Then the molecular operator in Eq. Q assumes the form 

d{t) = Y.K expHe,,(t-to)A] 



'eig-eo + i{TL + Tn)/2 



C£g(to) ■ (13) 



Upon substituting this result into Eq. ([S]), we obtain for 
the operators in the leads 

ctgit) = Q,(to)e-''^^'(*-*°)/'' 
+ 2^ 7— ^-T7rF-T-FrT7^^«'9'(*o) 



eo+i(rL+rR)/2 
X B[ei'q' -eeq] , 



i'q' ^ 



where, 



BiE) ^r{-]- inSiE) 



(14) 



(15) 
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FIG. 2: (color online) Currents (top row) and zero-frequency 
components of PSDs of accompanying noises (bottom row) 
for charge (left column) and heat (right column) transport 
through the single-orbital wire as functions of orbital energy 
for r = 0.1 meV. The remaing parameters are Tl = 5.2 K, 
Tr = 3.2 K (solid lines). For equal chemical potentials 
Ml ~ MR — the heat flow in panel (b) vanishes for equal 
temperatures Tl = Tr; its noise intensity in equilibrium at 
the temperature Tl = Tr = 4.2 K is depicted versus orbital 
energy £o by the dashed line in panel (d) . 



and V denotes the integral principal value [33| ■ 

Next we insert Eq. p^ and Eq. pi]) into the heat 
current operator, Eq. ([7]), and by consequently taking 
the ensemble average, we obtain a Landauer-like formula 
for the heat current [J, [23, [S^ ; reading, 

{J\t)) ■■= J'' = ^ J dEEr{E)[h{E)-ME)], (16) 

where the transmission coefficient, T{E) = Ti^r-n/KE — 
£o)^ + F^], is energy-dependent. Below we consider the 
case of symmetric coupling between the wire and the 
leads, Fl == Fr = F. 

The expression for the electric Scebeck current [J| reads 
very similar to Eq. ([T5)) , except for its absence of the en- 
ergy multiplier E in the integral in the rhs of Eq. (|16p . 
This seemingly small difference changes, however, the 
physics of the transport through the wire, because the 
multiplier inverts the symmetry of the integral. Namely, 
the Seebeck current is an antisymmetric function of or- 
bital energy and vanishes when the orbital energy level is 
aligned to the chemical potentials of the leads, Fig.[21[a), 
while the heat current is a symmetric function and ac- 
quires a nonzero value at Eq = 0, Fig.[2][b). 

The heat noise is described by the symmetrized auto- 
correlation function, i.e.. 



5'^(r) = l/2([Ajh(r),AJl'(0)]+), 



(17) 



of the heat current fluctuation operator AJ^(i) = Jg{t) — 
{jf{t)), where the anti-commutator [A, B]+ = AB + BA 
ensures the hermitian property. 

In the asymptotic limit t — > cx), the auto-correlation 
function depends on the time difference only. Its Fourier 



transform is the power spectrum for heat noise, obeying 
S^{oj) = S^{~uj) = / dTe^'^^S^ir) > , (18) 



A. The spectrum of heat fluctuations: explicit 

results 



Upon combining Eq. P^ and Eq. ([7]), a cumbersome 



being an even function in frequency and strictly semi- 
positive (Wiener-Khintchine theorem). In the following evaluation then yields the following, nontrivial explicit 
we address positive values of the frequency, lu > 0, only. expression for the PSD of electronic heat noise, reading: 

I 



s'\n = hij;n,Tj,) 




^±2 



O 



E±-] T{E)T{E±n) 



Tj [E{E - eo) -{E± n){E ± H - ep)]' 
[{E-eor + r^][{E±n~eo)^+T^] 



h{E)f^{E±n) 



E±^ T{E)T{E ± n)h{E)f^{E ±n) + [E^n{E)T{E ±n)T ^EnTiE)T{E ± n) 



n\ T'L'^T(E±n) - , 



n 



+ iE±n)[±j]T{E)r{E±n) 



-f 



TlT{E±n) 



^2"; (E-eo)^+T^ 



]ME)f^{E±n) 



r 



(19) 




FIG. 3: (color online) At zero-temperature Tl ~ Tr ~ the 
dependence of heat noise PSD versus frequency lj exhibits 
different power-law behaviors (dashed lines) for oj sampling 
intermediate values (proportional to uj^) as compared to the 
large frequency limit (proportional to a; ) when a; — >■ cxa. The 
other parameters are £o = and F = 0.1 meV. 



The PSD of heat noise at zero frequency uj — simpli- 
fies considerably, assuming an appealing form; reading 

5'^(c<. = 0;Tl,Tr) 

~JdEE^[TiE)[f^iE)[l-ME)] 

+ fn{E)[l - fn{E)]] 

+ riE)[l-TiE)][hiE)-MEr]. (20) 

This main result is in agreement with a conjectured pre- 
diction in Ref. [35|. The distinct difference between 
Eq. (PO]) and the PSD of the fluctuations of the nonlinear, 
accompanying Scebeck electric current, reading I2,l3|, 

5<='(c^ = 0;Tl,Tr) 
~ldE[TiE)[ME)[l-h{E)] 



+ /R(i?)[l-/R(i?)]] 

+ T{E)[l -T{E)][UE) - fn{E)]^] 



(21) 



wherein we abbreviated ^ = huj, f = 1 — f, and TZ{E) 
1 — T{E) is the reflection coefficient. 



is a factor E^ in the integral in Eq. ([20)1 . Although this 
distinction seemingly appears minor, it leads to a tan- 
gible difference in a dependence of the noise PSDs on 
energy level of the wire orbital, as depicted with Fig. [2] 
When the orbital level, Eq, is tuned to the chemical po- 
tential, the two expressions reveal different properties: 
while the zero-frequency component of the electric PSD 
at oj = exhibits a maximum at Eq = Oi its heat coun- 
terpart possesses a local minimum at this value. 



This difference originates from the saUcnt fact that the 
two transport mechanisms for charge and the energy are 
not equivalent. The electric current is quantized by the 
electron charge, e; in contrast, the energy carried by the 
electron is not quantized and may assume principally an 
arbitrary value. The main contribution to the electron 
flow across the wire stems from the electrons occupying 
energy levels around the chemical potential fi. 

Because the interaction between the leads and the 
molecule is weak, it only slightly perturbs the Fermi 
distributions, which possess strongly nonuniform profiles 
around /i = 0. Electrons of different energies contribute 
differently to the heat transport, but the Fermi distri- 
bution allows only for a finite number of electrons per 
energy level: i.e. just one in case of spinless electrons, 
or two in the case of electron with spin. Therefore, the 
electrons of given energy can move across the wire only 
when the corresponding level of the destination lead can 
host them. 

When £o deviates from the chemical potential, increas- 
ingly less electrons participate in the transport. The 
fiow of electrons becomes diminished, and since both, 
the electric current and the electric noise are insensitive 
to the electron energies, they both decrease with increas- 
ing 1^0 1- This scenario differs for heat transport: this is 
so because the deviation from the chemical potential in- 
creases the possibility that successive electrons will carry 
different energies. This in turn will lead to an increase 
of heat noise. With further deviation of the orbital en- 
ergy from the chemical potential, the occupancy differ- 
ence decreases monotonically and consequently the heat 
noise strength decreases again. 



B. Equilibrium heat fluctuations at Tl — Tr 

Next let us focus on the equilibrium properties of heat 
noise; i.e., the situation when the two temperatures are 
equal, Tl = Tr. In this case the average heat flow 
vanishes identically, but not its fluctuations. The zero- 
frequency spectra of both noises, for heat and electric 
noise, i.e. the corresponding power spectra increase with 
the increase of the coupling F, since it increases the trans- 
mission probability. The noise intensities are different 
from zero, however, even at equilibrium, see for heat noise 
Fig. [2j i.e., the two panels (b,d). 

The properties at vanishing temperature, Tl = Tr = 
0, are more subtle. It has been pointed in Ref. [S^ that 
heat noise seemingly violates the fluctuation-dissipation 
theorem (FDT). In fact, a temperature difference, AT, 
does not induce any thermal gradient. Therefore, there is 
no force which is conjugate to the heat current, and the 
process is out of the validity region of the FDT. This 'eva- 
sion' of the FDT is fully active in the zero-temperature 
limit, Tl = Tr = 0, where the heat noise PSD still de- 
pends on the frequency. This dependence is due to quan- 
tum fluctuations, the virtual transitions of electrons di- 
rectly from lead- to- lead [3^, [S^ . The Fermi distribution 



equals the Heaviside step function in this case. Thus, the 
contributions to the integrand in Eq. p^ comes from the 
interval [— ri, 0] only. After the integration of Eq. (|19p . 
one flnds the central result for the frequency dependent 
PSD: 



^h(c^,TL=TR 

r 



0) 



, [(217)2 -2r21 arctanf - 

4TTh r ^ V r 



2nT 



1 + log 



r^ 



{n^ + r^y 



, n = hLj. (22) 



In the limit F — >■ cxd the zero-temperature PSD thus scales 
like S^{u!) oc ui^. This is in full agreement with the results 
obtained in Refs. I3a,[33], where this asymptotic behav- 
ior has been found as being uniform throughout the whole 
frequency region. However, this is no longer the case 
when F is finite: the second term in the rhs of Eq. (j22[) in- 
troduces a linear cutoff in the limit w — s- 0, >S'^(w) oc uj. In 
distinct contrast, in the high-frequency region, the first 
term in the rhs of Eq. p2|) is dominating. As a result 
the PSD (P^ approaches a square-law asymptotic de- 
pendence, j{^) (X uj'^, in the high-frequrncy limit, see 
Fig. El 



IV. CONCLUSIONS 

By using the Green function formalism we have investi- 
gated electronic heat transport with our focus being the 
heat fiow fiuctuations for a setup composed of a single 
orbital molecular wire. For the noninter acting case we 
succeeded in deriving a closed form for the frequency de- 
pendence of heat current noise, i.e., the heat noise PSD, 
see Eq. (J19p . both in nonequilibrium Tl ^ Tr and in 
thermal equilibrium Tl = Tr. The dependence of the 
heat noise on the orbital energy Sq is qualitatively dif- 
ferent from that for the accompanying electric current 
noise. 

In the zero-temperature limit, the PSD of the heat 
noise obeys two distinctive asymptotic behaviors, being 
different in the low-frequency and the high-frequency re- 
gions. It is evident that the particular square-law shape 
of the PSD in the high-frequency region is due to the 
Lorcntzian shape of the transmission coefficient, T{E). 
Yet the general effect would remain for any choice of the 
coefficient in the form of a localized, bell-shaped function: 
the noise spectrum will deviate from a cubic power-law 
asymptotic behavior upon entering the high-frequency 
region. 

There is an intriguing perspective to apply an exter- 
nal periodic perturbation with the goal to control the 
spectral properties of the heat noise, similar to electronic 
shot-noise control in ac-driven nanoscale conductors [1^1 . 
This idea can be realized, for example, by subjecting the 
molecular wire to strong laser radiation [40| or by using 
direct modulations of the gate voltage. We conjecture 



that the role of laser radiation may give rise to novel 
phenomena to be explored further by combining a Flo- 
quet theory for the driven system with the nonequilib- 
rium Green function formalism |3|, llq . |33| . 

As emphasized in our introduction, with this work only 
the electron subsystem has been considered. Realistic 
heat transport in real molecular junctions would involve 
the complexity of interacting electrons and electron- 
phonon interactions [J] . As reasoned in the introduction, 
the electronic heat transport may dominate in certain sit- 
uations so that the measured heat noise can be attributed 
approximately to the electronic component only. The 
unified approach, which would include both the electron 
and the phonon subsystems, as well as the effects of their 



interactions, presents a future challenge, although several 
contributions in this direction for the heat current (but 
not the heat PSD) have already been undertaken before 
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